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Abstract 

In this paper, some new results are reported for the study of Riemann zeta 
function £(s) in the critical strip < Re(s) < 1, such as £(s) expressed in a 
generalized Euler product only involving prime numbers. Particularly, some new 
absolutely convergent series representations of £(s) based on binomial expansion 
are presented. The crucial progress is to find that £(s) can be expressed as a linear 
combination of polynomials of infinite degree, whose consequences are shown in 
several aspects: (i) numerically it provides a scenario to construct very fast 
convergent algorithm to calculate £(s); (ii) interestingly it shows that Lagrange 
interpolation using infinite number of integer Euler zeta functions reproduces the 
exact complex £(s); (hi) surprisingly it demonstrates that alternating Riemann 



zeta function (or other entire functions removing the pole of zeta function) is 
admissible to Melzak combinatorial transform for polynomials. Applying the 
functional symmetry on £(s) in the form of Melzak transform induces being 
written as the difference of two symmetrized factorials whose zeros are proved 
to all have real part of 1/2. Furthermore, the two symmetrized factorials are 
proved to have interlacing between the two sequences of the imaginary part of 
their zeros on upper (or lower) half plane, which ensures the difference of the two 
symmetrized factorials [proportional to ((s)] attaining the same feature of zeros 
with real part of 1/2 to endorse Riemann hypothesis. 
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1 Introduction 

Riemann zeta function ((s) is originally defined by analytic continuation via con- 
tour integral to extend Euler zeta function from the domain Re(s) > 1 to the whole 
complex domain^. If restricted on Re(s) > 0, ((s) has many equivalent representa- 
tions either by series summation or real variable integration. 

A popular series definition of ((s) on Re(s) > is via Dirichlet eta function rj(s): 



1 _ 1 y, (-l) n+1 _ 1 / _ 1_ 1_ 1_ \ 

n=l v / 



(1) 



where the series is conditionally convergent on Re(s)>0, and ((s) has a pole at s — 1 
[t)(s) is an entire function with 77(1) = In 2]. For s = 1 +i2mr/ In 2 (n 7^ 0) satisfying 
1 — 2 1 ~ s = were proved to be zeros of r](s) but not zeros 

of COO ^- All other zeros of 
r](s) are the same as those of C( s )- Riemann hypothesis is equivalent to state that any 
s on 0<i?e(s)<l which satisfies r](s) = must have Re(s) = 1/2 and Im(s)^Q. 
A well-known (Fourier) integral representation of ((s) on Re(s) > is 

C(s) = 7 ' ,x / dx = r- — r — — — / — e-^^dt (2) 

K J (1 - 2 1 -s)T(s) J e* + 1 (1 - 2^)T(s) J ^ + 1 K } 

There exist many integrals for ((s). Two examples are shown below (proofs omitted): 
fjn / r \ 

i n - 2+S 7T — T )dx=(-l) n (a-l)r(n-l+a)C(a) (neN, ite(a) > 0) (3) 







dx n \e x — l 



J x n - 2+s —\n^ Y —^jdx=(-irT(n-l + s)as) (neN , 0<Re(s) < 1) (4) 
We only consider Re(s) > due to the fact that ((s) has a functional symmetry: 



ai-s) = 2^^cos(^)T(s)as) (5) 



which can explore ((s) on Re(s) < by its symmetric counterpart on Re(s) > 0. 
Furthermore, on Re(s) > 1, Equation ([I]) becomes Euler zeta function which has no 
zeros proved by Euler product expansion. Thus only the critical strip < Re(s) < 1 
needs to focus where nontrivial zeros may exist. 

Theorem 1.1. A real positive s can not be a zero of Riemann zeta function such that 
C(s) ^ iflm(s) = and Re{s) > 0. 

Proof. From Eq. ^ on Re(s) > 0, if £(s) = 0, it must require that 

roo x Re(s)-l roa x Re(s)-l 

/ — cos[Im(s) \nx]dx = cos[/m(s) lnxi] / — dx = (6) 

J e x + 1 ' J e x + 1 

where mean value theorem for integral is applied, and x\ G (0, oo). As the latter integral 

is positive, it must have cos[im(s) lnxi] = that can not hold with Im(s) = 0. □ 



Theorem |1.1| is simple and known but very important result as it implies that 
any polynomial expression of ((s) will not have any real zero on Re(s) > 0. 

Theorem 1.2. The functional symmetry in Eq. |5p holds on 0<Re(s)<l. 

Proof. From another integral for ((s) on 0<i?e(s) <1 and series expansion of csch(x): 



we have the following identity: 

'x 



;i-2-)r( a )c( a ) = / ^(-^-i-^— d(|)=( 2 --i)C(i- a ) / ^ 

JO fc=1 ^ j +7r 2 J X +7T 



which is a simple proof of the functional symmetry in Eq. (|5J) on 0< Re(s) < 1. □ 
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Theorem L2 which is well-known since Riemann, is the key to answer Riemann 
hypothesis. The proof here is very simple, though only in critical strip. 



Theorem 1.3. ((s) on Re(s) > can be expressed into a generalized Euler product: 



\j=i 



Pi 



Pk + m X—B k 

i+hm J2 (aj-ir-^rK 1 -^ 

m— >oo ^— • X — S J - J - 

2j-X=p k+1 j=l 

2j-l\ 3,5,7,-, Ph 



(9) 



where {k, m} e N, and (2j — 1) includes the prime numbers larger than pk and all their 
possible products less than a sufficiently large prime number Pk+m- 



Proof. The summation-integral difference for n — > oo related to ((s) is 



111 1 

+ 3 1 + 5 1 + 7* + "' + (2ra- l) s A (2a;- 1) 



1 rfx = (l-2- s )C(s) + ' 



2(l-a) 



(10) 



Multiplying 3 s on Eq. (10) and subtracting from Eq. (10) yields 



(2n-l) s v 3Vi(2x-l 



: rfx=(l-3- s )(l-2- s )C( S ) + 



2(l-s) 



(11) 



Repeating the procedure till the kth prime number obtains a generalized Euler product: 



i + E - x ) 



i=(l+p fc+ l)/2 
2j-lt 3,5,7,-, p fc 



(2n - 1) 
1 - s 



l-s 



(12) 



3=1 



3=1 



An infinite large prime number Pk+ m replacing 2n — 1 turns Eq. (12) to Eq. (|9]). □ 

On Re(s) > 1, the summation and product in the square brackets in Eq. (JoJ) are 
towards zero when k — > oo which gives regular Euler product. On < Re(s) < 1, if k 
is chosen finite, terms in the square brackets in Eq. ^ converge after cancellation of 
divergent components; if k— >oo is chosen, and £(s)^0, the right-hand-side of Eq. ^ 



will be in the form of • oo in which L'Hopital's rule applies. If ((s) = £(1 — s) = 0, 



considered Eq. (12), the following equation holds for n— >oo: 

n 

1 + E w - !)" 



j=(l+p fc+ l)/2 
2j-lf 3,5,7,-, p k 



1 + 



1 - S 



(2n- 1) 



l-2_Re(s) 



(13) 



E - i) 5-1 

i=(i+p fc +i)/2 

2j-l f 3,5,7,-, p k 

which implies an equivalent Riemann hypothesis as: On < Re(s) < 1/2, if ((s) = 0, 



the norm ratio in the left-hand-side of Eq. (13) converges to 1 when n^-oo and k E N. 



2 Binomial Series Expansion of £(s) 

Convergent series expansions of ((s) are very limited. The following is such an 
example based on absolutely convergent binomial series: 

n=l v ' j=2 J \ J . 

Equation (14) is arrived by the series definition of i](s) in Eq. Q: 

1 oo 

V ( S ) = i + ^Et- 1 )^ 1 



n-1 



(Re(s) > 0) 



(14) 



i=2 



'( 8 










\d xj 







(15) 



*j=3 



followed by a binomial series expansion on (1 — X) s 1 with X = l— 1/j that converges 
absolutely for finite j, while may diverge for infinite j due to \X\ — > 1. Fortunately, due 



to the factor n/j arisen from the derivative, the series in Eq. (14) will converge 



Theorem 2.1. The binomial series expansion ofr](s) in Eq. (14) converges absolutely 
on Re(s)>0. 
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Proof. Comparison test of the series in Eq. (14) with Euler zeta function will be done 



as follows. For a sufficient large n, the absolute value of nth term in Eq. (14) is 



(n — s)\ 



n\(l-s)\ 



E 

3=2 



nj-iy- 1 / 1 

J 2 V J 



n—i 



o 



Re(s) 



n 



^ n(-l)H / l 
^ J 2 V J 



where the binomial coefficients asymptotic expansion 



n— 1 



(16) 



n — s 
n 



n s [-s)\ 



2/ n 



(17) 



is applied. For the summation of j in Eq. (16), we have 



E 

3=1 



n-l 



/7 



1 



(2j) 2 V 2j7 (2j + 1) 2 V 2j + 1 

oo , 

1 (2j) 2 J " U 



n-l 



<E 

3=1 



n 



(2j) 2 (2j + l) 2 



1- 



2j 



3 



n-l 



where the last step is realized by comparison with the integral result. Thus for suffi- 



ciently large n, considered Eqs. (18) and (16), it has 



(/>-*)! ^n(-l) j x ( l x 1 



TV. 



3=2 



r 



j 



o 



i 



l+i?e(s) 



n 



(19) 



which proves the absolute convergence of Eq. (14) on Re(s) > as compared to the nth 



term of Euler zeta function ((a) whose absolute convergence occurs on Re(a) > 1. □ 



In addition, many slightly different convergent binomial series expansion of ((s) 
or r](s) can be derived from a general form of 



;i-2 i - s )c( S 



s— a 

7 



s— a 

5 ~ I ■ 3=1 



E 



-IV'- 1 



/ 



Qq 

dx q j 



1- 1- 



3 / J 



(20) 



where q G N, {a, ^,7} G C with certain constraints for convergence. A special case of 
Eq. (20) on binomial series expansion of (1 — s)((s) has been studied in literature^. If 
choosing {a = 0, (3 — > 0, 7 = g = l}, it induces Taylor series expansion of r](s) at s — 1: 



■n(s) = 77(1) + lim ^ — 



1 \n / s—1 

1) / ^ 



n=2 v 7 j=2 



n-l 



n=2 v 7 *— n v / 



£=0 



m=l 



n=2 



(n-l)! 



-IV 



n-l 



(21) 



There also exist convergent binomial series of ((s) [instead of rj(s) in Eq. (14)] such as 



COO 



1 in ■ 

— + 

n=l v 



71 — S ! 



si 



1 n ( 1 1 

3=2 J V J 



n-l 



(ite(s) > 0) (22) 



But we will focus on the series expansion in Eq. (14), which is equivalent to 



00 / \ n 

X{n) 



n=2 ' j=2 

where x( n ) ( n > 1) is defined as 

00 , 



(ite(s) > 0) 



(23) 



i=2 



n-l 



(24) 



X{n) 



Some properties of x( n ) are x(l) = r /(2) — 1, = ^(1) — 1, and 



n=l 



1 



y(oo) = hm - Y^ 

3=1 



1 



' ji+i ' 



n-l 



n-l 







n 



(25) 



(^) (?) 
For a finite n, is the difference between midpoint and top-right corner rectangle 

approaches to the same integral, approximately exp(— l/x)/x 2 . As the sign of the 



difference summation in Eq. (25) depends on n, x{ n ) converges to zero in oscillation 
pattern [it is a subtle issue on the zeros of x( z ) as analytic continuation of xi n )}- 



Applying further binomial expansion on x( n ) m Eq. (24) will turn Eq. (14) into 

(26) 



v (n- 3 )! ^ (-l)^(fc + 2) 
V{ ) ^ (1 -s)\ 4^ k\(n-k- 1)! 



(Re(s) > 0) 



n=l 



fc=0 



where the binomial expansion on (1 -X)™- 1 with \X\ = 1/j < 1 and n e N absolutely 



converges, so does the series expansion in Eq. (26). It is worth to mention that numer- 



ically Eq. (14) is more favorable due to fast convergence of x{ n )i while Eq. (26) suffers 



a catastrophic cancellation problem in summation of k due to alternating binomial 
coefficients, which relies on high precision of inputs in order for accurate output. 

The crucial step towards revealing the zeros feature of £(s) is to switch the order 



of the two summations in Eq. (26) as follows: 



^ (n-s)! ^ (-l)N(fc + 2) 
n^L (1 - 8 )\ A? k\{n -k-l)\ 



n=i 



k=0 



lim ^ (-1)^ + 2) ^ 



[n — S ! 



k=0 



s)\k\ ^ (n-k-l)\ 

n=k+l V ' 



(27) 



which can be further written into 



lim ( m + 2 ~ s ) ! V ( m \ + 2 ) 

m^oo m\(l-s)\ ^\k) k + 2-s 

lim 7 ,,; TT j + 2 - s 

k=0 v 7 j'=0 

3+h 



(28) 



The validity of switching the two summations in Eq. (27) and the limit existence of 



m— > oo are ensured by absolute convergence of the two series. Equation (28) can be 



numerically verified. For example, a small m = 32 in Eq. (28) outputs a quite accurate 
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C(l/2) = -1.46034778. If {a = 0, (3 = 1, 7 = 2L, q = J} are chosen in Eq. pjh , a general 



version of Eq. (28) can be developed for rj(s) by using a subset of ((2k): 



7](s) = lim 



(m+l + J- 



s \ I "I 
2L 



-E 



2L^ fc=0 



m\ (-l) fc r?[2L(A;+J+l)] 



k + J+1- 



(JeN ,LeN) (29) 



2L 



2.1 ((s) calculated by fast convergent scheme. 



Equation (29) indicates that ((s) can be calculated by a subset of ((2k) with 
large k resulting in fast convergence. A numerical scheme to calculate ((s) with fast 
convergence can be constructed by using the following identity for ((s) with k E N, 
and Re(s) > 0: 

2k-l 



[l-(2k-iy-i}(T-l)((s) = J2 



n=2 



1 



(2*-i; 



n+l 



r(n + s ; 

2^n!f(V 



C(n + s) (30) 



Proof. Starting from the integral representation of ((s) 
"°° el - 1 



-x 5 " 1 ^ = (2 s -2)r(s)C(s) ( J Re(s)>0) 



(31) 



o 



we use geometric series summation and integral variable scaling in Eq. (31) to obtain 

2fc-l 



[i-(2fc-i) s - 1 ](2^-i)r( s )c( s ) 



2fc-l / j 



g 4fc-2 



J'=l 



e x ' - 1 



x s ^(fi (ite(s) > 0) (32) 



where fceN. Then Eq. (30) is arrived by taking power series expansion of exponential 



functions on numerator in Eq. (32) and integrating each term of the power series. □ 



In Eq. (30), k = 2 corresponds to the numerical scheme for calculation of ((s) as 



;i-3^)(i-2-K( s ) = i-| + ^ + ^-| + ^ 
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(33) 



which is arrived by multiplying 3 s on Eq. (10) and subtracting from Eq. (10). The 



series in Eq. (33) can be written into the form of Eq. (|30|) by binomial expansion: 
1 2 



£ 

k=l 



+ 



1 



(6A;-5) S (6A;-3) S (6fc-l) 



n=2 



T(n + s) (1-2 ■ 3" + 5 r - 



r(s) n\ 



-C(n+s) 



(34) 



where 1 — 2 • 3 n + 5 n = for n = 0, 1 leads to the convergence order of 0(n 1 2iie ( s )' 



roughly as fast as C(2)- In Eq. (30), k — 3 corresponds to 



14 111 1 4 1 

(1 - 5 1 - s )(l - 2- s )((s) = 1 + -- - + - + - + — + - _ + _ + ... 

K A ;sv ; 3 s 5 s 7 s 9 s 11 s 13 s 15 s 17 s 



(35) 



Linear combination of Eqs. (33), (35) and so on can generate series to calculate ((s) 
converging as fast as ({N), provided that proper combination coefficients (solved from 
a set of linear equations) let all terms of n < N vanish in the summation of n when 



combining k in Eq. (30). For example, the following combined series with the specific 



coefficients converges as fast as £(6) on Re(s) > 0: 



1 4 

l +n — ^ + 

3 s 5 s 



2 1 

3 I + 5 I + ' 


\ 10 s + 2 


) 756 


(3 s - 1 - 1) 


625 i 
"189 15 



18 s+1 / 1 1 1 8 
70 V 3 s 5 s 7 s 9 s 



81 
35' 



■vs-l 



(2 s - i)r(*)C00 



(36) 



It shows that Eqs. @, @, and Q need total 10 6 , 10 4 , and 85 terms [ 15, 25,45 terms 



chosen for the alternating series in Eq. (36) from left to right, respectively] to achieve 
0.37471336 - 0.27518432i, 0.36010325 - 0.26624621i, and 0.36010259 - 0.26624619, 
respectively, for calculation of C(0.2 + 2z) whose value is 0.36010259 - 0.26624620i. 
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2.2 rj(s) reproduced by Lagrange interpolation on a set of in- 
finite number of integer eta functions. 



Explicitly Eq. ( 28 ) can be reformed into 



which is an exact Lagrange interpolation formula on infinite number of integers 2, 3, 4, ■ ■ ■ 



The interpolation also can be done on a subset of integers based on Eq. (29). There 
exist infinite number of good table of nodes for convergent interpolation of r)(s). And 
some fast convergent iteration methods^ developed for Lagrange interpolation type 
equation can be applied to numerically find all the roots of r](s) = simultaneously. 

2.3 rj(s) is admissible to Melzak transform for polynomials. 

Melzak transform is inherited from combinatorial identities and finite difference 
theory for polynomials of finite degree. A basic Melzak transform is defined as: 

Theorem 2.2. If f(x) is a polynomial of degree m, the following transform holds 

/<*-^<V)g<Ml)T^r <^0,l,-,.») (38) 

for x, y e C Choosing x = yields a special case: 

/(-)=»( m ;")E(-i)*(:)s <^o,i,.., m ) (3 9) 

The proof is in literature—. Equation (28 ) shows that r](s) is admissible to Melzak 



transform by replacing y = s — 2, /(— k) —r){k + 2), and f(—y) = v( s ) m Eq. (39). Like 
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analytic continuation, Melzak transform can be used to extend a function from integer 
domain into complex domain (e.g., defining complex index Bernoulli numbers). 

Admissible condition to Melzak transform for m — > oo can be used to characterize 
an entire function which behaves like a pseudo-finite degree polynomial. If f m {—k) 
is the truncated polynomial of f(—k) by cutting off all degrees above m, then the 
requirement for f(—k) being admissible to Melzak transform in Eq. (39) for m— t-oo is 

"■n If ( m ~ ") B-D' (?) ~ / r ( " t)1 =0 (v * 0, 1, ... , m) (40) 

m->oo \ m J \ A; y y — k 

For instance, exponential decay and sine damping are admissible, which satisfy both 
lim f(—k) = and lim / (m) (0)/m = 0, while cosine and sine are not admissible. 

fe— >oo m— s>oo 

3 ((s) expanded by zeros of symmetrized factorials 



From Eq. (28), r](s) = is equivalent to 



E( ,,A (-1)^+2) 



m— ¥oo 

k=0 



k J k + 2 - s 

which unfortunately can not be solved analytically. To unveil the feature of r](s) = 



on 0<Re(s) < 1, we turn to apply various factorization by zeros on r/(s) in Eq. (28). 



3.1 r](s) as the summation of symmetrized factorials of all even 
degrees whose zeros all have real part of 1/2. 

The functional symmetry of ((s) or r](s) is the core reason to cause the special 

feature of ((s) = 0. We found that symmetrized factorials arisen from the functional 
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symmetry have all zeros with real part of 1/2. 

Theorem 3.1. All roots of the following symmetrized factorial polynomial equation 

n n 

TJ( a ._l + s ) + TJ( a ._ s) = (n>l) (42) 

have real part of 1/2 and nonzero imaginary part if all real (a-, — 1/2) do not vanish 
simultaneously and the nonzero (a,j— 1/2) have the same sign. 



Proof. Change the variable to be x = s — 1/2, Equation (42) becomes 

1 



n 

3=1 



a ,; V x 

3 2 



n 



1 

aj - - - x 



If any root x has nonzero real part as x = S + ir, then Eq. (|43|) requires 

i 2 n 



n 

3=1 



a,j - - + 5 + ir 
2 



n 

3=1 



a 7 — - + <5 I + r 



n 

3=1 



a,-, 5 — ir 

3 2 



n 

3=1 



(43) 



(44) 



However, when all (aj — 1/2) have the same sign, if S ^ 0, the norm ratio in Eq. (44) 



will always be greater or smaller than 1. Equation (44) holds only if 5 = or x = ir 



Thus all 2[n/2\ roots of Eq. (42) have real part of 1/2 and nonzero imaginary part 



Eq. (43) can not hold for real s = l/2 (i.e., x = 0) when at least one (aj — 1/2) ^ 0]. □ 



The anti-symmetrized version of Eq. (42) (the difference of the two factorials) 
can also be proved to have all roots with real part of 1/2 (and s = 1/2 is a root too). 
In this paper we only focus on the symmetrized version. 

Applying the functional symmetry of Eq. ^ on rj(s) in Eq. (23), we have 



n=2 



nl 



n-2 



n-2 



Y[(j + 2-s) + l[(j + l + s) 

.3=0 3=0 
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(45) 



which is valid in the critical strip because both r](s) and T](l — s) expanded via alter- 
nating series in Eq. Q are valid on < Re(s) < 1. Considered Theorem 3.1, Equation 



(45) can be factorized by zeros of symmetrized factorials of all even degrees: 



T](s) + T](l - s) 

3x(2) 



C(2) + 



2! 



£([i-(-i)i + "fV(-iri) ^ n ( 

n=3 V " ' 3=1 L V 



"J 



(46) 



Ad- 



here {1/2 ± iv n ,j} are the complex conjugated zeros of the symmetrized factorial of 



degree 2[(^ — l)/2j in the square brackets in Eq. (45). The first few {v n ,j} are 



{ I / 3j }=±^v / l5, {z/ 4j }=±^v / 7, Ki}=±^\A° 3± (47) 



with j = 1, 2, - ■ ■ , 2[(n — 1)/2J . If r)(s) is expanded by Taylor series in Eq. (21 ), since 
the zeros of (1 — s) n + s n are solvable, the factorization by zeros becomes 



rj(s) + 77(1 - s) 



LfJ 



n=2 



fe=l 



IV 1 

S -2 + i 



tan 



(2fc - 1)tt 
2n 



(4* 



A trivial case is Taylor expansion at s = 1/2 as factorization by a single repeated zero. 



3.2 r)(s) as the summation of symmetrized factorials of the 
same infinite degree whose zeros all have real part of 1/2. 



If the functional symmetry is applied on rj(s) in Eq. (28), it has 



T](s) + 77(1 — s) = lim y 



-l) k r](k + 2) 



«>- 00 - k\(m — k)\ 



j=0 



3=0 



(49) 
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Considered Theorem 3.1| r](s) in Eq. (49) can be factorized by zeros of the symmetrized 
factorials of infinite degree as 



■q(s) + 77(1 — s) = lim 



2(-l) k r](k + 2) 



n 



m->oo * — ' k\(m — k)\ 

k=0 v ' j=l 



(50) 



where {1/2 ± iT k j} are the complex conjugated zeros of each symmetrized factorial in 



the square brackets of Eq. (49), and an even m is chosen for convenience in this paper. 



3.3 r](s) as the difference of two symmetrized factorials of the 
same infinite degree whose zeros all have real part of 1/2. 



Equation (49) is a linear combination of the symmetrized factorials of the same 



degree. To deal with the case of linear combination, we have the following result: 



Theorem 3.2. The equation below is a linear combination of symmetrized factorials: 



in 

k=0 



Y[( aj -l + s) + l[(a 3 - s) 



j=0 



j=0 



(m > 1) 



(51) 



with all coefficients Ck of the same sign. Then all roots of Eq. (51) must have real part 



of 1/2 and nonzero imaginary part if all real (a 3 - — 1/2) do not vanish simultaneously, 
and the nonzero (a 3 - — 1/2) have the same sign. 



Proof. Equation (51) requires 



m m 



$^c fe JJ( Gi - s) 



k=0 j=0 



Ck 



3=0 



k=0 



a k - s 



m m 



(52) 



k=0 j=0 



3=0 



k=0 
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where the norm ratio should be 1: 



3=0 



E 

k=0 



Ck 



0-k ~ s 



H{ aj -i+s) 

3=0 



E 

k=0 



Ck 



a k - 1 + s 



(53) 



Assume that there exists a root s = 1/2 + 5 + it for Eq. (51). We can define the norms 



Vj = \dj - s\ 



Rj = \a,j — 1 + s\ 



(54) 
(55) 



for j = 0, 1, • • • , m. Then Eq. (|53|) becomes 

-l 2 




E 

.k=0 



Ck 
r 2 I 



1 



-5 



+ r^ 



-I 2^ 



E 

,fc=0 



Ck 




fc=0 
m— 1 m 



fc=0 j=0 



fc=0 £=fc+l 




n 

3=0 
3+k,+i 



m m 



m—1 in 



E^IW + EE 2 ^ 



fc=0 j=0 
3+k 



k=o e=k+i 



(56) 



n * 



Suppose that all real (oj — 1/2) do not vanish simultaneously, and the nonzero (a, — 1/2) 
have the same sign. If 5 has the same sign as the nonzero (aj — 1/2), it has rj < Rj for all 



j, then Eq. (56) can not hold for all Ck of the same sign because the numerator is always 



smaller than the denominator. If 5 and the nonzero (aj — 1/2) have opposite sign, then 
the numerator is always greater than the denominator. Thus if 3 = 1/2 + 5 + is a 



root of Eq. (51), it must have 5 = 0. And for real s = l/2 (i.e., r = 0), Equation (52) 

can not hold when at least one (aj — 1/2) ^ 0. □ 
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From Theorem 3.2 we have the following factorization by roots of Eq. (51): 

m 

/ m \ 2 

= 2 e<* n 



in 

k=0 



Y[( aj - s ) + Y[( aj -i + s) 



3=0 



3=0 



>fc=0 



s — 



T 



(57) 



where {1/2 + iTj} are the corresponding roots. On the other hand, from Theorem 3.1 
we have another factorization by zeros of each symmetrized factorial: 



m 
k=0 



Y[( aj -s) + Y[(a 3 -l + s) 



j=0 
3+h 



3=0 



E*hn 

j=l L 



k=0 



+ T k,j 



(5? 



where {1/2 + iTkj} are the zeros of the kth symmetrized factorial. Comparing Eq. (58) 



to Eq. (57) for s = 1/2 reveals one correlation among the imaginary parts of all zeros: 



k=0 j=l 
m 

k=0 



T, 



k.j 



(59) 



Specifically, in Eq. (49), if the same sign coefficients are grouped separately, The- 



orem 



3.2 can be applied on Eq. (49) to factorize by zeros of two combined polynomials: 



r/(s) + r](l — s) = lim 2 I \ 

m — ^rxn 1 < J 



r]{2k + 2) 



n 



t _ (2fc)!(m-2fc)l / _ 



-51 +e; 



E,, t k + 3 \ v)fl \( --)'+*! 

^ (2jfe + l)!(m-2Jfc-l)! I y [V 2J J 



(60) 



where {1/2 ± z0j} and {1/2 ± are the zeros of the combined polynomials of 



even and odd k terms, respectively. Equation (60) is arrived by first applying the 



functional symmetry on Eq. (28), and then combining the symmetrized factorials with 



coefficients of the same sign in Eq. (49). The order can be switched. Before applying 



the functional symmetry, we can first combine the factorials with coefficients of the 
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same sign in Eq. (28), or alternatively add an auxiliary function to convert all the 



combination coefficients in Eq. (28) to be positive and then subtract it [which avoids 
to deal with even and odd k terms]. For example, choosing f(x) = ( ^o^ 2771 ) / (^m™) ^ or 



Melzak transform in Eq. (|38j), then for m > 1, y ^ 0, 1, • • • , m, and x = 0, we have 

71—1 j 

( 61 ) 



m — y\ sr~- I 2k\ (2m — 2k\ 1 



E 



m / 14 — ' \ k 

k=0 



m — k J y — k 



2m\ (2m — 2y 
m ) V 2m y _ 2 2 " 
ml 



m — y 
m 



Substituting y = s — 2 in Eq. (61) reduces to 



7T 



m— 1 



[(' 



TT/ 5 N ^ 2- 2m 7r /2/c\ /2m — 2k\ 4-r. 

i^no+-,-*)=Y, W Ei W [ k )i m _ k )no+ 2 - s ) 



3=0 



fc=0 



(62) 



3=0 



which is a special case of the following identity from Melzak transform: 



7T 



m—1 



sin(7r/3) 



n^ +i -M=E l fe!(m %, no +7 



3=0 



fc=0 



(63) 



3=0 



where /3 £ Z and 7} G C. 



For < k < m, the summation coefficients in Eq. (62) are always positive: 



-2m „ 



71 (2k\ (2m - 2k 



[(^)!] 2 )\m-k 



1 (jfe-l)!(m-Jfe-i)! 



[( 



m-l \|]2 



fc!(m — k)\ 



> 



k\(m — k)\ 



(64) 



where we considered the fact that for < k < m, (™) has peak value at (^3) • 



Taking m— »oo in Eq. (62) and adding with Eq. (28) obtains 



ri(s) = lim > - 

m. — \r>c\ f J 



x . „ (A;-i)!(m-A;-f 
•l)^(A; + 2) + - ^ 



m—1 _ 
7T -n- , 5 



fc=0 



fc!(m — k)\ 



-no+»-)-T7^nw+5- 



lim $> fc TT(j + 2- S ) 



3=0 

m—1 K 
7T -n- , . 5 



3=0 



k=0 j=Q 

3+k 



[(¥)']' 3=0 
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(65) 



where all > due to Eq. (64) and r](k + 2) < 1. The sum of the coefficients is 



s, , , (k-\)\(m-k-\)\ 
■l) k r l (k + 2)+ y - 21 y 



[(^W X(m+1) , vr 



k=0 



k=0 



k\(m — k)\ 



(m+1)! [(m=l)\y 



(66) 



where x( m + 1) is defined in Eq. (24). 



Equation (65) can be numerically verified. For example, choosing a small m = 20 



in Eq. (65) outputs 0.70082616+0.43532002i for r/(0.2+2z) whose value is 0.70077353+ 



0. 43513124?. For combination with all €fc>0 in Eq. (65), we have the following result: 



Theorem 3.3. Assume that f(x) is a polynomial of degree m with positive leading 
coefficient, and have all real distinct zeros of {ai,02, ••• ,a m }. Polynomial g{x) of 
degree m—1 is obtained by linear combination as 



f(x) f(x) 
g(x) = ci^^ + • • • + c m ^^- 



(67) 



Then g(x) has m — 1 real zeros {61,62, • • • , 6 m _i} which interlace with m real zeros of 
f(x) as ai < 61 < 02 < ■ ■ ■ < a m _i < 6 m _i < a m if and only if all Ci are positive. 



The proof can be found elsewhere^. Applying Theorem 3.3 on Eq. (65) yields 
X(m+1) vr 



lim 

m— >oo 



m—1 



m—1 



(m+1)! I7m-i\n 2 



n(j +2+ dj -s) 5 TT (j - 
, f/ m-l Ni] 2 11 VJ 

LI 2 )■] J j=o IA 2 )■] 3=0 



= lim [f d (s) - f h (s)} 



(68) 



where < <i 3 < 1 is required by roots interlacing from Theorem 3.3 



Theorem |1.1| denies any real zero for r\ (s) on Re(s) > 0. Thus there can not 

exist any interlacing segment for three or more zeros between the zeros {2 + <i ,3 + 

c?i,--- ,m + l + d m -i} of fd(x) and the zeros {2.5,3.5, ■■■ , m + 1.5} of fh(s), otherwise 
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the combined polynomial will have at least one real zero. For example, if only one zero 
j + 2 + dj of fd(s) is in-between two consecutive zeros {j + 2.5, j + 3.5} of fh(s) to form 
a three- zero interlacing, then at two ends where fh{s) = 0, the combined polynomial 
fd(s) — fh( s ) will have opposite sign, the same as fd(s). This means that at least one 
time sign change occurs in-between (causing a real zero) for the combined polynomial. 
Considered both < dj < 1 and nonexistence of zeros interlacing segments, the only 
possibility to arrange the two sequences of the zeros of fd{s) and fh(s) to ensure the 
combined polynomial rj(s) having no real zeros on real positive s will be 

2.5 < 2+d < 3 + di < 3.5 < 4.5 < • • ■ < m+rf m __ 2 < m+l+d m _i < m+1.5 (69) 

where the smallest zero is from fh{s), since it has fh(s) + 1/2 < fd(s) < fh(s) + 1 due 



to 1/2 < r](s) < 1 for all real s > 0. Equation (69) also can be expressed as 



I TTl 

< d 2l+l <-<d n <l (£ = 0, 1, • • • , -) (70) 



We numerically verified the validity of Eqs. (69) or (70). Moreover, here we list a few 



properties of the two functions fd(s) and fh(x) defined in Eq. (68): 

(i) Two consecutive zeros of fd(s): {2£+2+d 2 g, 2£+3+d 2 £+i} with £g [0, m/2] are 
within the interval of two consecutive zeros of fh(s): {2£+2.5, 2£+3.5} with £g [0, m/2]. 

(ii) On real s>0, fd(s) and fh(s) cross the real axis m times, and /<f(s) is always 
on top of fh{s) more than 1/2 but less than 1 without crossing each other. 

(iii) fh(s) and fd(s) are monotonic decrease on s G [0, 2.5), and only have one peak 

between two consecutive zeros [as fjp(s) and (s) change sign one time in-between]. 

The following graph illustrates the displacement of fd(s) and fh(s) on real s>0: 
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Figure 1: Schematic graph of factorials fd(s) and fh(s) in Eq. (68) on real s > 0. 



Considering the functional symmetry of Eq. (68), we will have 

"to— 1 



7](s) + 77(1 — s) — lim 



I *( m+1 ) 1 



7T 



7T 



2 

TO— 1 



m— 1 



U=o 



i=o 

TO— 1 

UCjH-2-hrf,- — s)+ n^ + 1 + d i + s 

j=Q j=0 



(71) 



(m+1)! [(2^1)!]' 

which shows that r](s) [proportional to r](s)+r](l — s)} can be expressed as the difference 
of two symmetrized factorials of infinite degree that all zeros of each symmetrized 



factorial have real part of 1/2. In details, applying Theorem 3.1 on Eq. (71) obtains 



77 (s) + r/(l — s) = lim 



2tt 



2 X (m+l) | 2tt 



(m+1)! ' [(^)!] 2 ; n 



1 



-2) + ^ 



(72) 



In Eq. (72), ±Xj (j = 1, 2, • ■ ■ , m/2) are m real roots of the following equation 



TO— 1 



m— 1 



F fc (z) = l[(j + 2-ix) + n(./+2 + 



i=o 



3=0 



'TO— 1 



2 cos I arctan 



.r 



J + 2 



TO— 1 



n v / (j+2)2+x 2 =0 



(73) 
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which is arrived by substituting s = 1/2+ix (x is a real number!) in Eq. (71 ) as the zeros 



of the symmetrized factorials based on Theorem 3.1 Similarly ±u>j ( j = 1, 2, • • • , m/2) 



in Eq. ( 72 ) are m real roots of the following equation: 



m— 1 g m— 1 g 

j=0 j=0 



' m— 1 



2 cos | arctan 



a; 



j + | + dj 



m— 1 

n 



j + - + | x- 



(74) 



Below we will analyze the relationship between the roots of Eqs. (73) and (74). 



3.4 j]{s) as a single polynomial of infinite degree whose zeros 
all have real part of 1/2. 



Since Eqs. ( 74 ) and ( 73 ) are polynomials of only even power of x [all odd power 
of ix in Fh(x) and Fd(x) vanish otherwise Fh(x) and Fd(x) can not be real], we can 
define the new variable as y = x 2 so that {oo 2 } with j = 1, 2, ■ ■ • , m/2 are the roots of 
F d (y) = and {A^} with j = 1, 2, ■ ■ ■ , m/2 are the roots of Fh(y) = 0. Then we have 



Theorem 3.4. The m/2 distinct real roots {X 2 } of Fh(y) in Eq. ( 73) (defining y = x 2 ) 



interlace the m/2 distinct real roots {oj 2 } of F d (y) in Eq. ( 74) (defining y = x 2 ) or vice 



versa as X 2 < uj < X 2 < ■ ■ ■ < X 2 m/2 < u 2 m/2 or u 2 < X 2 < uj 2 < ■ ■ ■ < u 2 m/2 < X 2 m/2 



depending on specific sequence of0< d 2 e+i < 1/2 < d 2 e < 1 in Eq. (74) 



Proof. First of all, {X 2 } and {oj 2 } are all real guaranteed by applying Theorem 



3.1 



on 



Eq. (71). Finding the roots of F h (y) = [or F d (y) = 0] is equivalent to finding the 



roots of Fh(x) = [or Fh(x) = 0] on x > subject to a square mapping. 
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Since Fa(x) and Fh(x) have the same degree and the same sign of the leading 
coefficients, if Fd(x) and Fh(x) have sign interlacing on x > 0, then their zeros are 
interlacing on x > 0. Without loss of generality, we consider if there is a sign change of 
Fd(x) between any two consecutive real positive zeros of Fh(x). Let the two consecutive 



zeros of F h (x) on x > be x and xi, from Eq. (73) we have 



cos 



cos 



/m_1 / x \\ m ~ 1 ( x \ 2A; + 1 

> arctan ( 11=0 =>- > arctan I I = n (75) 

' m ~ x / j. \\ m ~ 1 / x \ 2k + 3 

> arctan I 11=0 ==>- > arctan ( ) = 7r (76) 



where fc G No. The sums of the angles in Eqs. (75) and (76) must have a difference of 
7r because from x to x\ is a continuous process and between x and X\, F^ix) (equiv- 
alently the cosine function) does not change sign (otherwise, they are not consecutive 
zeros) so that xq and X\ corresponds to two consecutive cosine zeros. 



If Fd(x) has sign change on x and x\, from Eq. (74) it must satisfy 

f m— 1 / \ \ /m—l 



COS 



arctan ( ^ | I cos f arctan ( ^ I I < (77) 



Considered the inverse tangent formula of sum of angles: 

m— 1 / \ m—l / \ m—l 



III, — X / \ III — -L / \ I III — X I \ 

> arctan ( 5 I = > arctan ( I — > arctan I rr (7£ 



with Pj(x) = (j + 2) + x , Equation (77) is equivalent to require 



■m—l / \ \ / m—l 



Sill 



arctan ( — -, — r- I I sin ( arctan ( — -, — r I I > (79) 



where Eqs. (75) and (76) are utilized. Equation (79) can be algebraically proved to hold 



for real positive x$ and x\ and dj satisfying Eq. (70). But it is more straightforward 
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to prove Eq. ( 79 ) by geometric graph showing that both sums of angles in Eq. ( 79 ) are 



less than tt/2. The following graph shows the complex vectors aj + ix: 



Im 



V2+d A 



5/2+d, 



7/2+d 2 
9/2+d 3 



AAA \ ; \ [A . ..A 




Figure 2: . Geometric graph of the angles in Eq. (79). 



in which the angles (containing the balls) are ZA OAi = arctan 

ZA 2 OA 3 = - arctan ( — X ° — — 

V3+pi(xo)/(di - 1/2) 



Xq 



2+p (x )/(d - 1/2) J' 
etc. It is obvious that even the sum of ab- 



solute angles in Eq. (79) are less than tt/2 so that Eq. (79) and then Eq. (77) hold, 
which proves the sign interlacing between F^(x) and Fh(x) on x > 0. Equivalently, it 
concludes that the zeros of Fh(y) and Fd(y) are interlacing, but the smallest zero can 
not be determined by a general sequence of < d 2 i+i < 1/2 < d 2 e < 1 unless additional 
information applies. □ 



In Eq. (72), when factorization by zeros between the two polynomials is done 



before the limit m — > oo is taken, then for a sufficiently large m, xi m + 1) can De 
chosen as a tiny negative number for convenience, which is similar to Taylor expansion 

of e~ x being terminated at a sufficiently large term with coefficient of either sign while 
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approaching the same convergent limit. If x(m + 1) < is chosen, then applying 



Theorem 3.4 will lead to an unambiguous zeros interlacing as Xf < co{ < \% < ■ ■ ■ < 



Kn/2 < u m/2 due ^° the requirement of 7/(1/2) > with x( m + 1) < in Eq. (72) 



Therefore, changing the variable as y — (s— 1/2) , Equation (72) indicates that 



T)(s) + r/(l — s) can be expressed as the difference of two polynomials of y 



7/(s)+7/(l — s) = lim 



(m+l) 



2n 2 )ff(y+^) - — ^ t\(y+ x l) 



JO) 



with interlacing zeros being — w^ /2 < — A^ 2 < • • • < — A 2 < — oj\ < — Af as a result of 



Theorem 3.4 and x( m + 1) < 0. 



For general linear combination of two polynomials, the following result is known: 

Theorem 3.5. Suppose that {a\, • • • , a n } and {bi, b 2 , ■ ■ ■ , b n } are all real zeros of 
f(x) and g{x), respectively, and g[x) interlaces f(x) as {b% < a\ < b 2 < • • • < b n < a n }. 
For a combined polynomial F(x) such that F(x) = af(x) + (3g(x) where a, (3 are two 
real numbers, if F(x) and g(x) have the same degree and have leading coefficients of 
the same sign, then F(x) has all real zeros {ci, c 2 , ■ ■ ■ , c n } and f(x) interlaces F(x) as 
{ai < Ci < a 2 < ■ ■ ■ < a n < c n }, provided that (3 < 0. 



A more general version was proved elswhere^. Equation (80 ) can be written into 
[r](s)+T)(l — s)]= lim 



2tt 



2 X (m+l) 



2tt 



m+l)! 



= lim 

m— >oo 



2tt 



G d (y) 



1) 



where the two polynomials Gh(y) and G^(y) have all real zeros {— A^} and {— wj}, 

respectively, and their zeros are proved to interlace as — cu^/2 < — < • • • < ~A 2 < 
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-UJ 



\ < —Xl- Moreover, the combined polynomial — [t](s) + — s)} and Gd{y) have 



leading coefficients of the same sign as the chosen x{ m + 1) < 0. In addition, the 



combination coefficient of Gd(y) is a negative constant. Therefore Eq. (81) satisfies all 



conditions to apply Theorem |3.5| to yield 

■m 

r)(s)+r){l-s) = hm - — 7Y^\\\y+^i) = hm 

m— >oo (m+lj! - LJ - rn— >ot 



(m+1)! 



n 



(82) 



3=1 v " ' ' 3=1 

where the zeros {— Qj} of t](s) + r/(l — s) are all real and interlacing with the zeros 



{-A^} of G h (y) in Eq. (81) as 



Am < -£lm < ■ ■ < -Xl < -0 2 < -A? < -Oj 



J3) 



Since {AH are all real positive numbers, Equation (83) indicates that {fij} also must 



be real positive numbers except for Q\ whose sign can not be determined by Eq. (83). 



If all Qj > 0, due to the negative leading coefficient [the chosen x(m + 1) < 0] in 



Equation (82), it will result in a contradictory 7/(1/2) < 0. Thus it concludes that all 



Qj > for j = 2, 3, • • • , m/2 except that Oi < 0. 

However, on the other hand, Q\ < means that r/(s) + r/(l — s) will have two 
real zeros as s = 1 /2 ± \/ — Oi, which seems to be contradictory to Theorem |l.l| In 
order to avoid this dilemma, it must have Q± < —1/4 such that the two corresponding 
"real zeros" s < and s > 1 are just out of the convergent domain (the critical strip) 
ofr/(s)+?7(l — s) as Eq. (82) is derived based on Dirichlet series in Eq. Q. Therefore, 



the factor of [(s — 1/2) 2 + Q x ] still exists, and Theorem 1.1 is not violated [within the 
convergent domain of r/(s) + r/(l — s)]. 
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Thus we have the final factorization by zeros for r](s) + T](1 — s) on < Re(s) < 1: 



T](s) + 7/(1 - s) 



(2 1 - - 2)tt- cos (y ) r(a) + (1 - 2 1 -) C, ( s ) 



lim !)■ 11 



m->oo (m + ; 



»4) 



where all > except that flj. < —7, and {fij} are determined by Eq. (71) as the 



difference of two symmetrized factorials. Equation (84) shows that ((s) is proportional 
to a single product of infinite number of quadratic forms [(s — 1/2) 2 + Qj] with all 
Qj > except that Q\ < —4, which immediately endorses Riemann hypothesis in the 
critical strip. 

It is worth to mention that 



(i) If x( m + 1) > is chosen in Eq. (72), it will also show that r](s) +T](1 — s) can 
be expressed into a single polynomial whose all roots have real part of 1/2. First of 
all, following the same procedure will arrive at the same conclusion that the zeros {A 2 } 
interlace the zeros {w 2 } or vice versa. In the case of oj\ < A 2 < uj\ < ■ ■ ■ < w^/ 2 < A 2 ^, 



Theorem 3.5 can be directly applied on Eq. (80) to conclude that the zeros {— uA 



interlace the zeros of {— Qj} leading to 

2x(m+l 



77(5) + 77(1 - s) 



lim 

m^oo (m + 1)! 



— TT 



3=1 



, - - j + 



with all {%} > 0. In the case of A 2 < u 2 < A 2 < ■ ■ ■ < \ 2 m/2 < 



m/2- 



55) 



defining the 



variable as y' = —(s— 1/2) 2 , Equation (72) becomes 



T](s)+T](l — s) 



- lim 

m— ¥00 



2 X (m + l) + 2tt 



2tt 



(m + 1) 



[(=?!)!]■ 



nw-*o-i7£biiw-«o 



3=1 



56) 
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Then Theorem 3.5 can be applied on Eq. (86) to conclude that the zeros interlace 
the zeros {flj} leading to 



, 2x(m+l) yj 2x{m+l) yj 

T]{s)+r](l-s)= hm M (0 3 - -y)= hm — — — M 

m-Vrm (^m+lj! - LJ - 



(m+1)! 



s ~2 J 



57) 



3=1 v ' i=i 

with all {fij} > 0. Thus regardless of the sign of x( m + 1)) v( s ) + vO- — s ) an d its 
approximation with finite large m will all have the expression of a single polynomial 
whose zeros all have real part of 1/2. 

(ii) Similar to Arndt-Gosper formula^, Fh(x) and Fd(x) in Eqs. (73) and (74) 
can be written into explicit forms of polynomial of x: 



F h (x) 



F d {x) 



2(-l) : 



m 



2(-l) : 



m 



TO— 1 TO— 1 

e n 

k=0 j=0 

TO— 1 TO— 1 

e n 

?c=o i=o 



., -(,; + 2) tan ( ^ 







3 , . / 7r(? — A;) 7r \ 
^ - (J + o + tan ( ^ + ^ ) 



m 



(89) 



for an even m. When m — > oo, under midpoint approach, Equations (88) and (89) 
become Cauchy principal integrals: 



F h (x) = lim 2(-l)T p.y. / [ f 



1 TO— 1 



j=0 
,1 TO— 1 

lim 2(-l)- P.!". / JJ 

i=o 



a; — (j +2) tan ( ——Try 



o 



^-(j + ^ + G^tan (~ _7r ^ 



= (90) 
rfy = (91) 



where all odd powers of x vanish due to 



P.V. 



, 71 J 
tan ( Try 

m J 



) dy 



(92) 



The interlacing relationship between the zeros {X 2 } of F h (y = x 2 ) and {u> 2 } of F d (y 



x ) may also be explored starting from Eqs. (88) and (89) 
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(iii) Equation (84) indicates that all roots of £(s) are a subset (equal to or less 
than the number) of the total roots {s = 1/2 ±i-y/Ttj}, since the prefactor of £(s) [i.e., 
k(s) = (2 1 ~ s — 2)7r~ s cos (tvs/2) r(s) + (l — 2 1_s )] might have isolated zeros. Nevertheless, 



Equation (84) implies that no matter where the roots comes from £(s) or the prefactor 
k(s), in the critical strip all their roots must have real part of 1/2. 

(iv) On the other hand, £(s) might have additional zeros (the trivial ones actually) 



that are not included in Eq. (84), because the expansion of £(s) in Eq. (84) is derived 
from Dirichlet series T)(s) whose valid domain is Re(s) > 0. Moreover, to make the 
functional symmetry 77(1 — s) also valid via Dirichlet series representation, the valid 
domain is restricted in the critical strip < Re(s) < 1. Thus all trivial zeros of £(s) 



located at real negative integers do not emerge in Eq. (84), and possible zeros of the 



prefactor k(s) which are outside the critical strip also will not show up in Eq. (84). 



4 Conclusions 

In this paper, based on absolutely convergent binomial expansion, alternating 
Riemann zeta function 77 (s) is found to be admissible to Melzak transform for infinite 
degree polynomials. Specifically, r](s) can be expressed as a linear combination of cyclic 



polynomials Pk(s) = J~| (j+2— s) with k — 0, 1, • • • ,m, which is shown in Eq. (28). 

3=0, 3+k 



Considered the functional symmetry of Riemann zeta function, the combined 
7](s) +77(1 — s), which is proportional to 77(5), can be written into a linear combination 



of symmetrized factorial polynomials P&(s) + -P&(1 — s), shown in Eq. (49). All roots 
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of Pfe(s) + Pfe(l — s) for each have real part of 1/2. Moreover, we proved that for 
a linear combination of Pfe(s) + Pfe(l — s) with same sign combination coefficients, all 
roots of the combined polynomial will still have real part of 1/2. 

Riemann hypothesis would be endorsed immediately if the linear combination 



coefficients in Eq. (49) all had the same sign. However, the combination coefficients in 



Eq. (49) are alternating between positive and negative, which equivalently leads to 7](s) 
expressed into the difference between two symmetrized factorials whose roots all have 
have real part of 1/2. Fortunately, we proved that the imaginary parts of the zeros 
of the two symmetrized factorials on upper half plane are interlacing. Based on well- 
known results about zeros feature of the combined polynomial from the difference of two 
interlacing polynomials, Riemann hypothesis is endorsed to show that the combined 
polynomial [proportional to C( s )] can De expressed into a single product of infinite 
number of quadratic forms (s — 1/2) 2 + Qj with all Qj > 0. 
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